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O ; ABSTRACT 



We report CCD photometric monitoring of the nonemission ultracool dwarf 



CLf BRI 0021-0214 (M9.5) obtained during 10 nights in 1995 November and 4 nights 

in 1996 August, with CCD cameras at 1 m class telescopes on the observatories 
of the Canary Islands. We present differential photometry of BRI 0021—0214, 
and we report significant variability in the J-band light curve obtained in 1995. A 
periodogram analysis finds a strong peak at a period of 0.84 day. This modulation 
^ ■ appears to be transient because it is present in the 1995 data but not in the 1996 

data. We also find a possible period of 0.20 day, which appears to be present in 
both the 1995 and 1996 datasets. However, we do not find any periodicity close 
to the rotation period expected from the spectroscopic rotational broadening 
(<0.14 day). BRI 0021—0214 is a very inactive object, with extremely low levels 
of H a and X-ray emission. Thus, it is unlikely that magnetically induced cool 
spots can account for the photometric variability. The photometric variability 
of BRI 0021—0214 could be explained by the presence of an active meteorology 
that leads to inhomogeneous clouds on the surface. The lack of photometric 
modulation at the expected rotational period suggests that the pattern of surface 
features may be more complicated than previously anticipated. 



Subject headings: — stars: low-mass, brown dwarfs — stars: rotation — stars: 
fundamental parameters — individual: BRI 0021—0214 
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1. Introduction 

During the last five years, we have witnessed the discovery of hundreds of ultracool 
dwarfs, 1 located in the solar neighborhood and in young open clusters and associations, e.g., 
Jones & Steele (2001). Detailed study of their properties is still in its infancy, but some 
general facts seem to be already well established. 

Non-gray model atmospheres have shown that dust forms under the low-temperature 
and high-density conditions prevalent in M dwarfs with T e ff < 2, 800 K, e.g., Tsuji, Ohnaka, 
& Aoki (1996a) and Allard et al. (1997). Dust grains such as calcium titanate (CaTi0 3 ) 
and corundum (A1 2 3 ) have been held responsible for the observed depletion of gaseous TiO 
from the atmospheres of the coolest M dwarfs. A new spectral class, designated as L, has 
been defined to classify dwarfs cooler than M (Kirkpatrick et al. 1999; Martin et al. 1999). 
The optical spectra of late L dwarfs do not show any TiO features, and they are dominated 
by broad resonance lines of alkali elements (Cs, K, Na, Rb). The J — K colors of late M and 
L dwarfs are extremely red due to their cool temperatures and the greenhouse effect caused 
by dust opacity, which dissociates the main molecular absorbers in the i^-band (CO and 
H 2 0). Dusty atmospheric models are more successful than dust-free models in reproducing 
the observed colors of late M and L dwarfs in the field (Leggett, Allard, & Hauschildt 1998), 
in the Pleiades cluster (Martin et al. 2000), and in very low-mass multiple systems (White et 
al. 1999; Leinert et al. 2000). On the other hand, the same models cannot explain the blue 
J — K colors of even cooler dwarfs (T e g < 1, 200 K), such as Gl 229 B, without gravitational 
settling of dust grains (Tsuji et al. 1996b). The real atmospheres of ultracool dwarfs are 
likely to be very dynamical, due to turbulent motions related to convection and rotation 
(Allard et al. 2001; Marley 2001). 

Ultracool dwarfs are generally fast rotators. Rotational broadening (v suu) for about 
40 late-M and L dwarfs has been measured using high-resolution spectra (Basri & Marcy 
1995; Martin et al. 1997; Tinney & Reid 1998; Basri et al. 2000; Basri 2001). The majority 
of the ultracool dwarfs have v suu > 5 km s _1 , with a trend for faster rotation toward cooler 
temperatures. The late L dwarfs also tend to be very inactive as judged from the flux of 
H Q emission (Gizis et al. 2000). Chandra observations of the M9 dwarf LP944-20 have set 
an upper limit to the quiescent X-ray emission that is two orders of magnitude lower than 
typical X-ray emission levels in M0-M6 dwarfs (Rutledge et al. 2000). 

BRI 0021—0214 is an ultracool dwarf (M9.5) that was discovered during the course 
of a survey for high-redshift quasars because it has an extremely red B — R color (Irwin, 
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McMahon, & Reid 1991). Among the known very late M and L dwarfs, BRI 0021-0214 is 
one of the fastest rotators with vsini = 37.5 ± 2.5 km s _1 (Basri 2001). Moreover, it is a 
very inactive dwarf because it usually does not show a detectable H a emission line (Basri 
& Marcy 1995), although a weak flare has been observed once (Reid et al. 1999). Even at 
the maximum of the flare, the level of activity was 3 times lower than the mean quiescent 
level in mid-type M dwarfs. The flare recurrence rate was estimated to be <7% of the time. 
Neuhauser et al. (1999) reported a nondetection of X-rays from BRI 0021—0214 in a pointed 
ROSAT/HRI observation of 106 kiloseconds, implying \ogL x /L ho i < —4.68, which is one of 
their most sensitive upper limits to X-ray emission from brown dwarf candidates. 

The fast rotation of BRI 0021—0214, originally reported by Basri & Marcy (1995), 
prompted us to start a program of CCD monitoring in 1995. Earlier results of a program of 
CCD observations of M7-M9 field dwarfs have been reported in Martin, Zapatero Osorio, & 
Rebolo (1996). Thermal inhomogeneities in the atmosphere, such as magnetic cool spots or 
massive cloud systems, could be modulated with the rotation period, which is expected to 
be shorter than 0.135 day for a typical radius of 0.1 R & . Thus during one observing night, 
the object completes a few rotation cycles. This paper reports on the results of monitoring 
the behavior of BRI 0021—0214 in the J-band over 10 nights in 1995 and 4 nights in 1996. 
In §2 we describe the data acquisition, reduction, and analysis techniques. In §3 we present 
the light curves that result from the differential photometry. We discuss the results and 
implications in §4, and our conclusions and directions for future research in §5. 



2. Observations and Data Analysis 

CCD images of BRI 0021—0214 were obtained on the nights of 1995 November 8-15 
and 22-23 with the 1 m Jacobus Kapteyn Telescope (JKT) at the Observatorio del Roque 
de los Muchachos (La Palma), and on the nights of 1996 August 19-22 with the 0.82 m 
telescope (IAC80) at the Observatorio del Teide (Tenerife). We used standard broadband 
/ filters. The JKT was equipped with a 1024 x 1024 pixel Tektronix camera mounted on 
the Cassegrain focus, providing a field of view of 32.3 arcmin 2 (0 / .'333 pix -1 ). The IAC80 
was equipped with a 1024 x 1024 pixel Thomson camera mounted on the Cassegrain focus, 
providing a field of view of 54.4 arcmin 2 (0 / /433 pix -1 ). All nights were photometric with the 
exception of 1995 November 15 and 1996 August 21-22, which had some patchy thin cirrus. 
Exposure times ranged from 300 s to 500 s. We adjusted the exposure time to compensate 
for extinction due to high airmass or cirrus clouds. The seeing of the images varied from l'.'O 
to 2 / .'5, and the airmasses at which they were taken ranged from 1.18 up to 1.80. All the CCD 
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frames were bias subtracted and flat fielded using packages within the IRAF environment. 2 

We performed aperture photometry on BRI 0021—0214 and on five other stars in the field 
using "vaphot," an IRAF script written by Hans Deeg (Deeg et al. 1998). We verified that the 
two brightest stars in the field of view are not variable on the timescale of our observations, 
and we have used them as reference stars. They are shown in Figure 1. Images were first 
aligned, and the nonoverlapping parts were removed. Changes in the seeing conditions and 
in the telescope focus throughout a night implied variations in the size of the point-spread 
function (PSF) of the stellar sources from frame to frame. To take this into account we 
derived an average PSF of the three stars by fitting a circular Gaussian on each frame. All 
three sources comply with the requisites for a good definition of the PSF. Additionally, to 
obtain high signal-to-noise aperture photometry, we determined optimum apertures for each 
star according to Howell (1989). All apertures were expressed in multiples of the FWHM of 
the average PSF, and the sky intensity was defined as an outer ring 1.5 pixels wide. These 
multiples were kept constant from night to night. Finally, instrumental magnitudes were 
computed, and relative magnitudes (i.e., differential photometry) of BRI 0021—0214 with 
respect to the sum of the flux of the two reference stars were derived. The 1 a accuracy 
of the differential photometry is ±0.006 mag and ±0.003 mag for the IAC80 and JKT data, 
respectively. 



3. Results 

Figure 2 displays the light curve of BRI 0021—0214 obtained at the JKT in 1995 Novem- 
ber. The photometry of reference star 1 with respect to reference star 2 is also shown. The 
target presents a larger photometric dispersion than the two reference stars. In Figure 3, 
we show the standard deviation (a) associated with the total JKT differential light curve of 
BRI 0021—0214 and five other stars (including the references 1 and 2) in the field of view of 
our CCD images. We call this quantity 07 in analogy to Bailer- Jones & Mundt (2001). The 
/-band photometry was calibrated using the DENIS J-band magnitude of reference star 1 (X. 
Delfosse, 2000, private communication). The value of 07 increases toward fainter magnitude, 
as expected. We find a rather tight relationship between /-band magnitude and 07 value 
for the JKT data (07 = 0.17144 - 0.028001 x I + 0.001154 x I 2 ) with a scatter of 0.0006 
mag (1 a). The amplitude of photometric variability in the JKT data of BRI 0021—0214 is 
0.0182 magnitudes. This value is 14 a away from the relation between I-magnitude versus 



2 IRAF is distributed by National Optical Astronomy Observatory, which is operated by the Association 
of Universities for Research in Astronomy, Inc., under contract with the National Science Foundation. 
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a i for the reference stars. 

Figure 4 displays the total light curve of BRI 0021-0214 obtained at the IAC80 in 1996 
August. The scatter in the differential photometry of BRI 0021-0214 in the IAC80 data set 
is a i = 0.006 mag, which is a factor of three lower than that in the JKT data set (07 = 0.018 
mag). On the other hand, the photometric errors in the IAC80 data are larger due to the 
smaller aperture of the telescope and poorer image quality. In Figure 5 we plot the 07 versus 
/-band magnitude (DENIS system). The relation between these two quantities is not as 
good as for the JKT data. The standard deviation of the fit shown in Figure 5 is 0.00316 
mag (compare with 0.0006 mag for the fit shown in Figure 3). This is probably because there 
are large PSF variations across the CCD frames taken with the IAC80. BRI0021-0214 and 
reference star 1 have a better PSF than the average of the other reference stars because they 
are closer to the center of the chip. Our differential photometry was obtained by averaging 
the fluxes of the reference stars, and this procedure gives more weight to the brighter sources. 
The reference star 1 is the brightest source. Figure 5 indicates that photometric variability 
of BRI 0021-0214 with the same amplitude as in the 1995 JKT dataset (07 = 0.018 mag) 
can be ruled out in the IAC80 dataset with a confidence of 2.5 a. Thus, the photometric 
variability of BRI 0021—0214 in the I-band appears to have decreased from 1995 November 
to 1996 August. In the next subsection we will show that this effect can be accounted for 
by a source of periodic modulation which was present in 1995 but vanished in 1996. 

Since BRI 0021—0214 is redder than the reference stars, we checked for the presence of 
correlations between differential magnitude and airmass or seeing, which, if present, might 
indicate color-dependent differential photometric effects due to refraction and scattering in 
the Earth's atmosphere. We did not find any such correlation, as illustrated in Figures 6 
and 7. We conclude that the variability of BRI 0021—0214 is likely to be intrinsic to the 
object, and not due to changes in the Earth's atmosphere. 

3.1. Periodogram analysis 

We performed a period analysis of the data using a CLEAN algorithm modified by one 
of us (H. J. Lehto, 2002, in preparation). For an earlier version of this method, see Roberts, 
Dreher, & Lehar (1987). The sampling in both the 1995 November and 1996 August data 
sets is such that the window function shows high sidelobes at small multiples of about 1 d _1 . 
In principle, these could create ambiguity in the determined period. 

Because CLEAN (and any other deconvolution method) contains inherently degenerate 
solutions, we ran several solutions with different values of the gain to get an idea of the stabil- 
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ity of the solution. The solutions for our periodograms remained stable. In 1995 November, 
BRI 0021—0214 shows a single strong peak with a frequency of 1.1898 d -1 , corresponding to 
a period of 0.8405 day. This peak has a pseudo-power about 2 orders of magnitude higher 
than the highest peak in the comparison star's CLEANed periodogram (Fig. 8). The ampli- 
tude of this peak is 0.018 mag, which is similar to the value of 07. In 1996 August there is no 
sign of the 0.8405 day peak, but a new small cluster of peaks has appeared around 4.8842 d _1 
(corresponding to a period of 0.205 day). These peaks clearly have a lower amplitude than 
previously, and may represent a sinusoid-like variation with a modulated amplitude. There 
is also a small possibility that the dominant frequency of this variation is actually 4.0115 d _1 
(corresponding to a period of 0.249 day). The amplitude of the peak at P = 0.205 day is 
0.007 mag, which is similar to the photometric errors. The periodogram can, however, find 
a sinusoid of amplitude below the photometric uncertainty because the noise is spread over 
many frequencies in the power spectrum. 

Figures 9 and 10 present the JKT light curve of BRI 0021—0214 folded with a period of 
0.84 day. In the first figure, data points taken at different nights are given different symbols. 
In the second one the data points have been binned in groups of five and a sinusoid fit is 
overplotted. We note that around phase=0.5 the data seems to deviate from the sinusoidal 
function. Each individual error bar corresponds to the standard deviation of a group of five 
points. Figure 11 displays the IAC80 light curve of BRI 0021-0214 folded with a period 
of 0.20 day. Data points have been binned by 5, as in Figure 10. The amplitudes of the 
sinusoidal fits are 0.022 mag for the JKT and 0.006 mag for the IAC80. After subtracting 
the 0.84 day sinusoidal fit from the JKT data, we searched for the 0.2 day period. The JKT 
residual photometry is plotted in phase with P = 0.2 day in Figure 12. Data points have 
been binned in groups of seven to increase the S/N ratio. The shape of the phased light curve 
suggests that the 0.2 day period may be present in the JKT data, although the amplitude is 
lower than the 0.84 day modulation. In fact, the amplitude of the 0.2 day sinusoidal fit to 
the JKT residuals is 0.006 mag, which coincides with the amplitude of the fit to the IAC80 
data. 



4. Discussion 

The presence of the 0.2 day photometric modulation in both the IAC80 and JKT light 
curves of BRI 0021—0214 suggests that this is a long-lived modulation. Could it be the 
rotation period? The absolute magnitude and T e g of BRI 0021—0214 are Mk = 10.15 and 
2,500-2,000 K (Leggett et al. 1998), respectively. Another observational constraint is that 
lithium has been depleted from the atmosphere (Basri & Marcy 1995). According to the 
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evolutionary models of Chabrier et al. (2000), the lack of lithium rules out the possibility that 
the object is younger than 0.3 Gyr and less massive than 0.06 M Q . Given these constraints, 
the models predict that the radius of BRI 0021-0214 should be R < O.lli? . For R = 
O.lli? , sini = 1, and V cq = 35 km s _1 , we get a period of 0.14 day, which is inconsistent 
with our periodogram analysis. The radius would have to be 0.14 R & to get a rotation period 
of P = 0.2 day, which is inconsistent with the theory. We consider several hypotheses to 
explain this discrepancy. 

First, we consider the possibility of non-rotational line broadening. Three different 
studies have derived wsini independently by using high-resolution optical spectra (Basri & 
Marcy 1995; Tinney & Reid 1998; Basri 2001). They all find values of about 40 km s _1 , and 
we have adopted the most recent one. The agreement between v sini measurements obtained 
at different epochs indicates that the broadening is not due to an unresolved spectroscopic 
binary. BRI 0021—0214 does not have an unusual rotation among the sample of late M 
and L dwarfs that have been studied at high spectral resolution. It is unlikely that the 
broadening is due to any line formation effect, which should basically be the same for all 
dwarfs of similar gravity and temperature. 

Second, we consider the possibility that the radius is about 30% larger than predicted 
by the models. Recently, the radius of an extrasolar giant planet has been observed for 
the first time. The planet HD209458b has a radius about 40% larger than that of Jupiter 
(Mazeh et al. 2000). The favored interpretation is that the radius has been puffed up by 
stellar insolation (Burrows et al. 2000). The radius of BRI 0021—0214, however, cannot be 
increased in the same way because of the lack of a close main-sequence star. The lack of 
detected X-ray emission by ROSAT (Neuhauser et al. 1999) rules out the possibility that 
there is a compact relativistic primary, such as a black hole or neutron star. The effects of 
rotation have not yet been included in models of brown dwarfs. One possible effect in a fully 
convective object is a larger radius and a cooler central temperature (Martin & Claret 1996). 
Preliminary calculations by I. Baraffe (2000, private communication) indicate, however, that 
the rotational distorsion in a fast-rotating 0.1 M star is small (<1%). 

Our previous discussion has ruled out the possibility of problems with either the vsini 
or the theoretical radius. Thus, we must address the evidence that BRI 0021—0214 presents 
periodic photometric variability with periods longer than the rotation period. 

Cool stars are usually understood by analogy to the Sun. Periodic photometric variabil- 
ity in cool stars arises from the presence of dark spots, where the enhanced magnetic field 
suppresses convective energy transport and cools off the plasma to temperatures below that 
of the photosphere. Early- and mid-M-type dwarfs have a high occurrence of H Q emission, 
light variability, and flares (the BY Dra syndrome), e.g., Bopp & Espenak (1977); Herbst & 
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Miller (1989). The photometric periods of BY Dra stars are consistent with their expected 
rotation periods. The fraction of the bolometric luminosity that is emitted in the H Q line 
(i?H Q = Ln a /L ho \) is log i?H Q = —3.9 ± 0.2 for the active, photometrically variable, BY Dra 
stars. Many of these stars have photometric variability with an amplitude of 0.3 to 0.02 mag 
in the V-band (Bopp & Espenak 1977; Pettersen et al. 1987; Alekseev & Gershberg 1997). 

BRI 0021—0214 is rotating even faster than most BY Dra stars, but it is not active. 
It has logi?H a < —6.3 (Basri & Marcy 1995), which is more than two orders of magnitude 
lower than the BY Dra stars. Magnetic flares cannot account for the observed light changes 
because they happen <7% of the time and are too weak to produce continuum light at 
optical wavelengths (Reid et al. 1999). Thus, it seems unlikely that the source of photometric 
variability in BRI 0021—0214 is a magnetic field because most of the time there is no evidence 
for chromospheric or coronal cooling. 

Bailer-Jones & Mundt (2001) recently presented results of an /-band search for vari- 
ability in 21 late-M and L dwarfs. Half of their objects present evidence for variability with 
amplitudes of 0.01 to 0.05 mag on timescales between 0.4 and 100 hr. CLEAN periodograms 
found strong periodicities at about a few hours in some of the objects, but for other objects 
the variability did not seem to be periodic. For one object with observations obtained 1 
year apart, no common period could be found. Bailer- Jones & Mundt (2001) argued that 
the periodicities were due to long-lived surface features modulated with the rotation peri- 
ods, and the absence of periodicities for other stars was due to surface features that evolved 
in brightness and physical size on timescales of a few hours to a few days. The surface 
inhomogeneities could plausibly be due to magnetically induced spots or plages, or dusty 
clouds. 

If magnetic fields produce surface features, we would expect that the amplitude of the 
photometric modulation may scale with the strength of H a emission. We could also expect 
that photometric modulation becomes smaller for cooler objects, because they tend to be 
more inactive (Gizis et al. 2000). Bailer- Jones & Mundt (2001) found that the amplitude of 
photometric modulation does not depend on spectral type from M6 through L5. In Figure 13 
we plot photometric variability (07) versus H a equivalent width (EW Hc J of BRI 0021—0214 
together with L dwarfs for which Bailer- Jones & Mundt found evidence for variability. For 
these ultracool dwarfs there is no correlation between EWh q and the amplitude of photo- 
metric variability, suggesting that these two observables are not causally connected. 

We conjecture several possibilities to explain the pattern of variability in BRI 0021—0214: 
(a) There is differential rotation in BRI 0021—0214, however it would have to be ex- 
treme to explain the very long periodicity of 0.84 day; (b) There are turbulent regions in 
BRI 0021—0214 reminiscent of Jupiter's bands and zones. Perhaps bright spots are lagging 
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behind the mean rotation of the photosphere; however, the period of 0.84 day seems again 
to be too extreme because Jupiter's zonal velocities are moderate (Smith & Hunt 1976). 
On the other hand, this period may be related to a zonal oscillation. Jupiter's Great Red 
Spot has a zonal oscillation with a period of three months (Solberg 1969), and it is possible 
that an ultracool dwarf may have shorter zonal oscillation period; (c) The fast rotation of 
BRI 0021—0214 (almost three times faster than Jupiter) may produce instabilities related to 
strong Coriolis forces. The instabilities may lead to wave motions typical of rotating fluids, 
similar to inertial waves or Rossby waves e.g., Tritton (1988). It is possible that the 0.2 day 
periodicity corresponds to a stable pulsation mode, but the 0.84 day periodicity cannot be 
accounted for by stable pulsations; (d) Convection reaches the optically thin layers of ultra- 
cool dwarf atmospheres (Hauschildt, Allard, & Baron 1999). Coriolis forces and convective 
overshooting may push clouds of dust some distance beyond the theoretical radius of struc- 
tural models. Some cloud decks could reach out to significant altitudes with respect to the 
surface. They would be dragged by rotation, but they could lag behind the deeper layers. 
The photosphere of BRI 0021—0214 may not rotate as a solid body. The 0.2 and 0.84 day 
periods could arise from bright or dark clouds located high above the photosphere. 



5. Conclusions and Future Directions 

We have monitored the brightness of the ultracool dwarf BRI 0021—0214 for 10 nights 
in 1995 November and 4 nights in 1996 August using /-band CCD imaging. BRI 0021—0214 
has shown significant variability, particularly in 1995. We find a surprisingly long modulation 
in the 1995 data set (0.84 day) that cannot be due to rotation. A CLEAN periodogram of 
the 1996 data set gives a periodicity of 0.2 day. None of these two periods correspond to the 
expected rotation period (<0.14 day). 

Since BRI 0021—0214 is a very inactive object with extremely low levels of H a and X- 
ray emission, we favor the interpretation that the source of surface thermal inhomogeneities 
required to explain the photometric variability are clouds rather than magnetically induced 
spots. BRI 0021—0214 has weak TiO bands, which is a sign of Ti depletion by condensation 
onto grains. The atmosphere of BRI 0021—0214 is too warm for water to condense, but 
silicate and iron clouds are expected (Chabrier et al. 2000; Ackerman & Marley 2001). 
Theoretical models indicate that they play a key role in controlling opacity and temperature 
structure in the atmosphere. The rapid rotation and cool temperature of BRI 0021—0214 
could make its surface appear more analogous to Jupiter than to the Sun. Gelino & Marley 
(2000) have simulated integrated light observations of Jupiter at 400 nm (5-band) made by 
the Hubble Space Telescope. They find that the apparent magnitude could be modulated 
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with Jupiter's rotation period. The amplitude of the sinusoidal light curve would be about 
0.04 mag, which is similar to the amplitudes found by us and Bailer- Jones & Mundt (2001). 
Gelino & Marley (2000) predict that the amplitude would be larger in the thermal emission 
regime. Thus, near-infrared monitoring of BRI 0021—0214 and other ultracool dwarfs may 
be worthwhile. 

Tinney & Tolley (1999) reported a search for variability using a tunable filter technique 
to measure changes in the strength of TiO bands. Tinney (2000) discusses more techniques, 
such as long-slit spectra and narrowband imaging. We suggest that fT-band observations 
could be particularly valuable to search for changes in the dust-induced greenhouse effect on 
the CO bands. Optical spectroscopic monitoring could test the hypothesis that the clouds 
are made of dust because the strengths of the TiO bands should be variable. A report of 
possible variability of the Li I resonance line in Kelu 1, an L2 dwarf, was made by Martin et 
al. (1998). Since the Li I line is formed against a background of TiO bands (Pavlenko 2000), 
the apparent variability of the atomic line could be due to variable TiO opacity and/or 
variable continuum opacity. Kirkpatrick et al. (2001) reported variability in the L dwarf 
Gl 584 C in the 835-975 nm region in timescales of days and years. 

We have briefly discussed four possible scenarios to explain the variability observed 
in BRI 0021—0214: differential rotation; bands and zones moving with different velocities; 
pulsations; and turbulent clouds located high above the photosphere. In any case, the 
variability of BRI 0021—0214 is more complicated than previously anticipated. Simultaneous 
multiwavelength observations may shed light on what phenomena are responsible for the 
surprising photometric behaviour of BRI 0021—0214. 
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I filter, 5.4' x 5.4' 




Fig. 1. — /-band image of BRI 0021—0214 obtained with the IAC80 telescope. The reference 
stars are indicated. 
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Fig. 2. — Complete JKT light curve of BRI 0021—0214. The lower panels display data 
obtained in individual nights. Open symbols are used for BRI 0021—0214 and small dots for 
reference star 1. 
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I magnitude 

Fig. 3. — Standard deviations versus /-band magnitudes for BRI 0021—0214 (filled circle) 
and 5 reference stars (open triangles) in the JKT database. A second-order polynomial fit 
to the a i values of the reference stars is shown with a solid line. BRI 0021—0214 has a 07 
value significantly higher than the solid line, clearly indicating variability. 
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Fig. 4.— Complete IAC80 light curve of BRI 0021-0214. The lower panels display data 
obtained in individual nights. Open symbols are used for BRI 0021—0214 and small dots for 
reference star 1. 
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I magnitude 

Fig. 5. — Standard deviations versus /-band magnitudes for BRI 0021—0214 (filled circle) 
and 5 reference stars (open triangles) for the IAC80 data. A second-order polynomial fit to 
the 07 values of the reference stars is shown with a solid line. 
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Fig. 6. — BRI 0021—0214 data as a function of airmass. No correlation is seen. 
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Fig. 7. — BRI 0021—0214 data as a function of seeing. No correlation is seen. 
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Fig. 8. — CLEANed periodograms of BRI 0021—0214 and comparison star 1. The abscissa 
in each figure in the inverse of period units of d -1 . Panel (a), at the top, is the CLEANed 
periodogram of BRI 0021-0214 during 1995 November (JKT). Panel (b) is the CLEANed 
periodogram of comparison star 1 during 1995 November (JKT). (Note the change in the 
ordinate representing the pseudo-power.) Panels (c) and (d) are the respective CLEANed 
periodograms for BRI 0021-0214 and star 1 during 1996 August (IAC80). Panels (e) and 
(f) show the window functions of the 1995 November and 1996 August observing sessions. 
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Fig. 9.— Phased JKT light curve of BRI 0021-0214 with P = 0.84 day. Different symbols 
denote data points obtained in different nights. 
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Fig. 10.— Phased JKT light curve of BRI 0021-0214 with P = 0.84 day. Data points have 
been binned in groups of five. 
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Fig. 11.— Phased IAC80 light curve of BRI 0021-0214 with P = 0.20 day. Data points 
have been binned in groups of five. 
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Fig. 12.— Phased JKT residual light curve of BRI 0021-0214 with P = 0.20 day after 
subtraction of the P = 0.84 day sinusoidal fit shown in Figure 10. Data points have been 
binned in groups of seven. 
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Fig. 13. — Amplitude of variability in the J-band of ultracool dwarfs studied in our work 
and in that of Bailer- Jones & Mundt (2001) versus EWn a - Open circles denote H a emission 
detections. Filled circles denote H Q emission upper limits. 



